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The integrated Sachs- Wolfe effect (ISW) can be an important factor in the generation of Cosmic Microwave 
Background anisotropics on all scales especially in a reionized curvature or lambda dominated universe. We 
present an analytic treatment of the ISW effect, which is analogous to thick last scattering surface techniques 
for the Doppler effect, that compares quite well with the full numerical calculations. The power spectrum 
of temperature fluctuations due to the small scale ISW effect has wave number dependence k~ 5 times that 
of the matter power spectrum. 
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I. INTRODUCTION 



It is well known that several effects contribute to Cosmic Microwave Background (CMB) fluctuations, 
e.g., the difference in the gravitational potential between the last scattering surface (LSS) and the present 
(the Sachs- Wolfe effect) [1], Doppler shifts from moving electrons on the LSS [2], and the fluctuations 
corresponding to adiabatic or isocurvature initial conditions. The gravitational redshift induced by a time 
varying potential, i.e., the integrated Sachs-Wolfe (ISW) effect, gives another important contribution to 
CMB fluctuations although it vanishes in linear theory for an fio = 1 flat universe. This effect is dominant 
on very large scales [3], but only the non-linear effect is usually considered important on small scales [4]. 
In this paper, however, it is shown that the ISW effect does indeed lead to significant small scale CMB 
fluctuations for curvature or cosmological constant dominated models. In particular, under the assumption 
of early reionization, this effect is quite important since the Doppler contributions are severely cancelled 
on scales below the thickness of the LSS [5]. Such early reionization is actually required for open primeval 
isocurvature baryon (PIB) models [6]. We introduce here a simple model-independent analytic method for 
calculating this effect which agrees quite well with the detailed numerical solution. 



II. ANALYTIC FORMALISM 

The full first order Boltzmann equation [7] for the evolution of temperature perturbations 0(x, 7, 77) in 
the gauge invariant notation [8] is given by 




where Vi is the matter velocity (c = 1), 11^ is the anisotropic stress perturbation, 7, are the direction cosines 
for the photon momentum, T is the Christoffel symbol, over-dots are derivatives with respect to conformal 
time rj = J dt /a, &o is the monopole component of 0, and \f is the Newtonian potential. The last term in 
eqn. (1) accounts for Compton scattering, where f = x e n e <TTa is the differential optical depth, with x e the 
ionization fraction, n e the electron number density, and a? the Thomson cross section. The perturbation 
to the intrinsic spatial curvature <J> is given by the Poisson equation (V 2 — 3A')<J? = 4irGa 2 pS, where the 
curvature constant K = —Hq(1 — fio — A), the Hubble constant Ho = lOO/i km s _1 Mpc -1 , the scaled 
cosmological constant A = A/3_0~q, and the total density fluctuation 6 = Sp/p. 

In the absence of scattering, the right hand side of (1) vanishes for a flat (flo = 1) universe since <J> = 
— = in first order. We then obtain the usual free streaming equations and the ordinary Sachs- Wolfe effect, 
i.e. 0(x,7, 770) + *[x(»7o)) Vo] = 0[x(r?*),7, 77*] + \F[x(r7*), 77*], where 770 and 77* are the present and the last 
scattering epoch respectively. The more conventional expression 1 /3\D r (fc, 77*) is obtained from a combination 
of and \f on the LSS (see [7]). Moreover, the quantity + \f is the temperature anisotropy apart from 
the unobservable present day monopole and dipole. In open and (flat) lambda dominated universes, \t — $ 
does not vanish and an additional effect, the ISW effect arises. 

It is possible to obtain a simple approximate analytic solution to eqn. (1) by analogy to thick-LSS 
techniques [9, 10]. In a pressureless matter dominated universe, <J> = — \t [7], and on small scales, we may 
drop the curvature terms. We can then solve this equation in Fourier space after Compton drag has become 
negligible z < z d ~ 130(fi /j 2 ) 1/5 ^ 2/5 : 

0(fc, p, 7? ) + 770) ~ [0(fc, p, Vd ) + Vd )}e ik ^-^e-^^ + Qssf + &isw, (2) 

t(?7i, 772) = T(rj)dr) is the optical depth, kp = k • 7, ®ssf is the secondary scattering fluctuation (SSF) 
generated on the new LSS if there is reionization, and &isw is the ISW effect. Explicitly, these terms are 
given by 

Qssf = I °(©o + * + pv)i e-^'^e^^-^dr], 

V (3) 

Q ISW =2 / ^(k,r])e- T (™o) e ii=»(v-vo) dT]j 
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where &ssf has contributions from the Doppler [9] and ordinary Sachs- Wolfe effect. 

Note that eqns. (2), (3), and the following argument are valid for any ionization history. In particular, 
for standard recombination x e (z < zj ~ 1000) = 0, and r(r? ( j ) %>) = t{t], r)o) = &ssf = 0, i.e. only the ISW 
effect adds to ©(fc, ji, r)d), the primary temperature fluctuations after standard last scattering. On the other 
hand, for reionized scenarios, the primary fluctuations are damped by exp(— r) due to rescattering [11]; in 
no-recombination scenarios, they may be ignored since r(r? ( j ) *7o) ^> 1- 

Whereas the integrand of &ssf peaks at the LSS, the integrand of &isw is localized near the curvature 
dominated epoch (1 + z g ) ~ l/fio — 1> or lambda dominated epoch (1 + z g ) ~ [1/fio — l] 1 ^ 3 . This epoch 
defines what may be called a gravitational "last scattering" surface (GLSS). A simple and useful analytic 
solution for &isw can be obtained on scales smaller than the thickness of the GLSS by noting that the LSS 
cancellation arguments for the SSF terms (valid on scales smaller than the thickness of the LSS) [5,11] can 
apply to the GLSS as well. 

Extending the techniques of Kaiser [9] and Efstathiou [10] to the ISW effect, we find that the present 
day power spectrum of temperature fluctuations 

k 3 IV P(k) I' 1 

k 3 P(k) = —|0 + ~ -^-ttH / \Gssf(x) + G ISW (x)\ 2 dx, (4) 

where |© + , i r |J! J7ls = 1 fe dyu|©(fc, r]o) + 1 ^(k, r]o)\ 2 , x = r]/r]o, the matter power spectrum P(k) = (\S(k)\ 2 ), 

Gssf(x) = \tiI(Dt + bT)e- r <™°\ 

3 ( 5 ) 
Gisw(x) = -a- 2 r,lHfo (ba - Dh)e~ T ^o\ 

a = 1/(1 + z), and D is the growth factor in linear theory [12] such that 6(r/) = D(r/)6(rio). 

Notice that the wave number dependence of the ISW contribution is identical to the secondary Doppler 
term: V(k) oc k~ 5 P(k). This dependence is easy to understand physically. Only the [i = mode in O 
can survive cancellation across a thick LSS or GLSS since the contributions are then additive. In general 
this restriction suppresses fluctuations by a factor (fcr?o) -1 ^ 2 - However, for the Doppler and ordinary Sachs- 
Wolfe contributions, there is an additional suppression by a factor (fcr?o) _1 since both sources in the original 
equation (1) are proportional to [i. Generalizing refs. [9, 10], we have set f(©o + + fiv) = ifc _1 [iif + 
vt+0(t 2 ^) + 0(f V?)]. Since \f oc v/k, the ordinary Sachs- Wolfe term is negligible on small scales compared 
with Doppler term. On the other hand, the integrated Sachs- Wolfe term, which does not vanish for the 
jj, = mode, has the same dependence in k as the "cancelled" Doppler term. It is therefore not negligible 
on small scales. 

III. COMPARISON WITH NUMERICAL RESULTS 

For comparison, we solve eqn. (1) numerically up to present by the method outlined in [13]. It is 
instructive to examine first the time evolution of the quantities in eqn. (1). In Figure 1, we have chosen 
a representative no-recombination model (see caption) to illustrate the evolution of \f, v, and |© + Vflrms 
in (a) an open and (b) a lambda (A = 1 — fio) model. After the drag epoch z < 100, v and \f evolve in 
linear theory. At the curvature or lambda dominated epoch, \f starts to decay, marking out the GLSS. More 
precisely, whereas the LSS is defined by the visibility function fexp(— r), the GLSS is given by \D r exp(— r) 
(dashed lines, arbitrary normalization). 
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FIG 1: Evolution of \f, v, and the r.m.s. of + \f with and without ISW from numerical calculations for 
a specific small scale mode (k = 0.03Mpc _1 ) in (a) an open model (A = 0) and (b) a lambda model 
(A = 1 — flo) f° r an = = 0-1) h = 0.5 universe with no recombination and isocurvature initial 
conditions. Overall normalization is arbitrary. Through the LSS, temperature fluctuations, originally 
of O(v), become increasingly suppressed due to cancellation and rescattering. We have also plotted the 
shape of the LSS [oc fexp(— r)] and the GLSS [oc \D r exp(— r)] from linear theory for reference (dotted 
lines). In lambda models, the ISW effect contributes strongly near the present epoch, and the GLSS is 
not well defined. Nevertheless the analytic approximations still work quite well. 
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FIG 2: Power spectra for the temperature fluctuations k 3 V(k) where k 2 = k 2 + K [3]. from analytic 
and numerical calculations: (a) No-recombination; open (upper set) and lambda (lower set) baryon 
isocurvature models with spectral index n = and parameters fio = &b = 0.1, h = 0.5 where the 
initial entropy fluctuations |5(fc)| 2 oc k n . The solid lines represent the numerical results with (thick) 
and without (thin) ISW. The dashed lines represent analytic calculations in an analogous manner, (b) 
Standard recombination; open (upper set) and lambda (lower set) CDM models with adiabatic initial 
conditions P(k) oc k n and parameters fio = 0.1, = 0.03, h = 0.5, n = 1. On large scales, the power 
spectra fall off since we have subtracted the unobservable monopole and dipole. 
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We have plotted the resultant temperature fluctuation |0 + \t|,. ra> with (thick line) and without (thin 
line) the ISW term in eqn. (1). Notice that the ISW has a negligible effect during the last scattering epoch; 
only in the lambda or curvature dominated epoch and after suppression through the LSS does the ISW effect 
become important. In these no-recombination scenarios, the LSS and GLSS have comparable width, and 
the ISW effect is important across the spectrum. 

We plot the full present day power spectrum of temperature fluctuations with a as = 1 normalization 
for various models in Figure 2. In panel (a), the results for no-recombination isocurvature open (top lines) 
and lambda models (bottom lines) of Figure 1 are shown. The analytic calculations (dashed line) with 
(thick) and without (thin) ISW agree quite well with the full numerical solution on small scales. Deviations 
on large scales are expected since the analytic formalism is only valid below the thickness of the LSS and 
GLSS. We calculate the ratios of the power spectrum of temperature fluctuations with and without the ISW 
contribution for other baryonic no-recombination models in Table 1. Notice that the ISW effect gives a 
contribution of the same order as the secondary Doppler effect for low fio models. 

TABLE I 
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TAB 1: Ratios of the power in temperature fluctuations with to without ISW R = l^tot /T^noi sw given 
by the numerical (R num ) and analytical (R an ai) computations for an = fio, h = 0.5 universe. 
Rnum is calculated for k = lfio^Mpc -1 but does not vary significantly on small scales. The ISW 
effect is less significant in lambda models since lambda domination occurs at a relatively late epoch. 
Analytic predictions agree quite well with numerical results (< 20% in power, < 10% in temperature 
fluctuations). 

In Fig. 2b, we examine the ISW effect in an open (top lines) and lambda model (bottom lines) with 
standard recombination. Here we take an adiabatic CDM model for which this ionization history is more 
natural. Although this model overpredicts fluctuations if normalized to <r%, we choose it to demonstrate that 
even in this extreme case, the ISW effect plays no role on small scales. For standard recombination, the LSS 
is relatively thin, and the ISW effect is masked until the Silk damping scale. Furthermore since the GLSS 
is much thicker than the LSS, cancellation of the ISW effect begins at a much larger scale and therefore is 
negligible at the Silk damping scale. Note however that the ISW effect itself does not depend sensitively 
on the ionization history. Moreover, even in these models, the ISW plays an important role on large scales 
(where the analytic formalism breaks down). It is interesting that the analytic formalism nevertheless 
predicts the correct scale for dominance of the ISW effect. Finally, we should mention that the ISW does 
contribute significantly to open or lambda adiabatic CDM models with sufficiently early reionization or no 
recombination. 
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IV. DISCUSSION 



We have found that the ISW effect is important at all scales in a reionized open or lambda universe 
including the so-called Doppler peak and below. Furthermore, we have developed an analytic formalism 
which compares quite well with detailed numerical solutions. From this formalism, we find that the power 
spectrum of temperature fluctuations behaves on small scales as V(k) oc k~ 5 P(k) for both the Doppler and 
ISW contributions, despite the fact that the scale dependence of v and \f themselves are not similar, due to 
a difference in their cancellation behavior. The amplitude of the temperature fluctuations can be predicted 
to 10% accuracy by the analytic approximation. In standard recombination models, the small scale ISW 
effect is masked since the LSS is relatively thin compared to the GLSS, and thus cancellation of the Doppler 
effect does not occur until a very small scale. The analytic approach derived in this Letter is a potentially 
powerful tool for understanding and calculating the ISW effect. It is independent of the detailed model and 
is thus applicable to all scenarios which have a time varying potential, e.g., even topological defect models 
[14] for which the ISW effect is essentially the only important contribution to CMB fluctuations. 
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TABLE CAPTIONS 

TAB 1: Ratios of the power in temperature fluctuations with to without ISW R = l^tot /T^noi sw given 
by the numerical (R num ) and analytical (R an ai) computations for an = fio, h = 0.5 universe. R num 
is calculated for k = lfio^ 2 Mpc _1 but does not vary significantly on small scales. The ISW effect is less 
significant in lambda models since lambda domination occurs at a relatively late epoch. Analytic predictions 
agree quite well with numerical results (< 20% in power, < 10% in temperature fluctuations). 

FIGURES: 

FIG 1: Evolution of \f, v, and the r.m.s. of + \f with and without ISW from numerical calculations 
for a specific small scale mode (k = 0.03Mpc _1 ) in (a) an open model (A = 0) and (b) a lambda model 
(A = 1 — flo) f° r an = = 0-1) h = 0.5 universe with no recombination and isocurvature initial 
conditions. Overall normalization is arbitrary. Through the LSS, temperature fluctuations, originally of 
O(v), become increasingly suppressed due to cancellation and rescattering. We have also plotted the shape 
of the LSS [oc fexp(— r)] and the GLSS [oc \D r exp(— r)] from linear theory for reference (dotted lines). In 
lambda models, the ISW effect contributes strongly near the present epoch, and the GLSS is not well defined. 
Nevertheless the analytic approximations still work quite well. 

FIG 2: Power spectra for the temperature fluctuations k 3 V(k) where k 2 = k 2 + K [3]. from analytic and 
numerical calculations: (a) No-recombination; open (upper set) and lambda (lower set) baryon isocurvature 
models with spectral index n = and parameters fio = &b = 0.1, h = 0.5 where the initial entropy 
fluctuations |5(fc)| 2 oc k n . The solid lines represent the numerical results with (thick) and without (thin) 
ISW. The dashed lines represent analytic calculations in an analogous manner, (b) Standard recombination; 
open (upper set) and lambda (lower set) CDM models with adiabatic initial conditions P(k) oc k n and 
parameters fio = 0.1, fi^ = 0.03, h = 0.5, n = 1. On large scales, the power spectra fall off since we have 
subtracted the unobservable monopole and dipole. 
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